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ABSTRACT 



Four modes of possible rotor motion are defined: synchronous preces- 

sion (without rub), partial rub, full annular synchronous rub and reverse 
whirl. A model is developed to investigate the relationships between a 
rotational system's characteristic parameters and the conditions of equili- 
brium for the full annular rub motion. The dynamics of an unconstrained 
rotor with unbalance are first established. The impact of constraining the 
rotor's motion by a rigid casing with a finite clearance is then analyzed. 
For any given system, the speed regimes in which full annular synchronous 
rub is in equilibrium are defined as a function of the system's unbalance, 
clearance, rotor stiffness, damping and the coefficient of friction between 
the rotor and casing. The rigid stator restriction is subsequently allevi- 
ated and the equilibrium conditions for full annular synchronous rub are 
reanalyzed. Additionally, reverse whirl as a possible stable mode of rotor 
motion is studied. Specific rotational systems are then analyzed according 
to the model developed to demonstrate the effect of varying particular 
parameters. 
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INTRODUCTION 



All rotating machinery have some practical degree of unbalance which 
drives flexural vibrations in the system. These vibrations may be readily 
calculated in terms of the system 1 s mass, unbalance, rotational speed, 
radial stiffness and damping imposed by the medium in which the rotating 
element operates. Should the magnitude of the vibrations exceed the 
physical constraints of the enclosure which houses the rotating element, the 
potentially destructive phenomena described here as rotor rub occurs. 

Shaft vibrations, misalignment, blade creep and thermal expansion are 
some of the physical realities which precipitate design of a radial 
clearance between the high speed rotor and surrounding casing in 
turbomachinery . Theoretical efficiencies of the turbomachinery are reduced 
though, by blade tip losses which result from the presence of the radial 
clearance between the rotor blade tip and the casing. Demand for improved 
system efficiency thus drives a reduction of the blade tip clearance and so 
increases the likelihood of rotor rub. Rotor rub in fact may be a design 
condition and so accounted in the design by the use of abradable materials. 
[ 1 ] 

Prior to the occurrence of rub, the rotating element’s dynamics may be 
described as synchronous precession, driven by the system’s residual 
unbalance. If rub develops, it may occur over a fraction of the rotor’s 
precessional orbit (termed partial rub), or may persist over the full orbit 
(full annular rub). Of these two conditions, full annular rub is 
potentially more destructive and is the subject of this thesis. For further 
information regarding partial rub, attention is drawn to reference [2-4]. 
Examples of full annular rotor rub occur with rotor blade-to-casing or 
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rotor-seal rub, as well as with a rotating shaft having clearance in a 
bearing. [5] The rub introduces a reaction force at the rotor’s periphery 
which strongly influences the rotor’s dynamics and may additionally excite 
the casing (further referred to as the stator) into motion. The reaction 
force may be broken down into its normal and tangential (friction) force 
components; the two being linearly related by the coefficient of friction. 
Friction may be moderated by proper lubrication (such as in the case of 
bearing rub) as well as through suitable selection of materials. Other 
design considerations or abnormal operating conditions may preclude these 
preventive measures though, and so friction may take on a particularly 
destructive nature. 

Friction is the resistive force which acts tangentially at the 
interface between two surfaces in relative motion. By acting tangentially 
on the rotor’s periphery, it essentially establishes a torque which counters 
the driving (input) torque for the rotating element. The result is either a 
loss in rotational speed or a higher input torque required to maintain a 
constant speed. The effect of friction on the mating surfaces is to cause a 
local heat buildup and, in the very least, some degree of grinding wear and 
deformation of the rotor and seal or bearing annulus. The more extreme case 
encompasses catastrophic destruction of the machine. [2] Typically 
accompanying these mechanisms is the annoying and often unacceptable 
acoustic phenomenon due to frictional vibration instabilities known as stick 
slip. Dynamically, the occurrence of rotor rub may destabilize the 
synchronous precession of the rotor. If friction forces are of large enough 
magnitude, the counter-torque developed may instigate a reverse whirl 
condition wherein the precession is in a direction opposite to the direction 
of rotation of the shaft. The nature of reverse whirl is such that the 
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rotating element is subjected to high frequency alternating stresses which 
may lead to rapid fatigue failure of the shaft. [6] 
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1 . DEVELOPMENT OF THE MODEL 



The characteristics of rotating machinery designed for the generation 
or transmission of power vary greatly with their particular application. An 
exact analysis of the motion is extremely complex, but many meaningful 
results can be determined through employment of reasonable simplifying 
assumptions in a single degree of freedom model. Even with these simplifi- 
cations though, numerous variables are necessary to describe the motion 
involved . 

The model in this thesis is displayed in Figure 1. The rotating 
element consists of a central disk of mass, M , mounted symmetrically on a 
laterally flexible but torsionally rigid, massless shaft of stiffness K . 

The shaft stiffness is determinable knowing the shaft’s material type, 
cross-sectional properties and length. The shaft is then supported at each 
end by rigid bearings without clearance. The natural frequency of the 
rotor, the frequency of excitation at which resonance occurs, is given 

by = /K /M^. The rotor disk is enclosed in its casing, separated 
radially by a symmetric clearance, r , existing when the rotor is at rest. 
(The effects of gravity are neglected in this analysis). The casing, or 
stator, is characterized by its mass, M^, and its mounting stiffness, K . 

As in the case of the rotor, the natural frequency of the stator is given by 
ui. = /K /M . Damping is introduced for motion of the rotor and stator 

INI S S S 

respectively by the damping coefficients, D and D . It is the presence of 

r s 

damping in both the rotor and stator which limits the magnitude of their 
responses at their respective resonant frequencies. Damping also causes a 
phase lag between the system's excitation and response. For ease of 
analysis, the model has been "compressed" in Figure 2 to a planar diagram 
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with the element’s dynamic characteristics represented as in a simple linear 
resonator system: a spring-mounted mass with a dashpot in which the energy 

of the motion is dissipated. Both radial and longitudinal symmetry are 
implicit in this model. 
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FIGURE L ROTOR RUB MODEL 
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FIGURE 2. RDTDR RUB PLANAR K2DCL 
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2. FREE ROTOR DYNAMICS 

It is assumed that under design conditions, the rotor is driven at a 
constant speed, to, independent of the external load applied. This assump- 
tion essentially neutralizes the effects of fluctuations in the design load 
due to damping or rub related friction. Under ideal conditions, it is 
imaginable that the rotor is perfectly balanced. Thus when driven at its 
design speed, the rotor center remains at the center of rotation; aligned 
with the axis defined by the supporting bearings' centers. This idealiza- 
tion is not realized in practice though, and to some extent, a residual 
unbalance exists in the mass distribution of the rotor. This unbalance may 
be depicted in the planar model by locating the mass center (or center of 
gravity), g, at a distance 6 from the center of the disk. The effect of the 
unbalance is to establish a centrifugal force acting upon the center of 
gravity while the rotor is in motion. This centrifugal force, which is 
proportional to the square of the rotational speed, acts as the driving 
force for the rotor's dynamical motion. 

An analysis of the motion of the free rotor (i.e., without displacement 
constraints imposed by the stator) provides the baseline condition for which 
the rub condition may be compared. The equation of motion for the free 
rotor is derived by applying Newton's 2nd Law to the rotor: 

(sum of forces acting on the rotor) = 

(mass of the rotor) x (acceleration of the rotor center of gravity). 
Referring to the free body diagram for the rotor, Figure 3. the forces which 
must balance in the steady motion of the rotor are described in the 



following: 
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FIGURE 3. FREE-RD70R FREE BODY DIAGRAM 
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a. Centrifugal Force - As described above, this force results from the 
presence of a residual unbalance; i.e., the center of gravity does not 
coincide with the center of rotation. The centrifugal force acts on the 
center of gravity and is directed radially outward from the bearing axis 
(stationary reference position about which the rotor precesses). 

b. Spring Force - The spring (or restoring) force acting radially 
inward on the rotor center results when the rotor center is displaced from 
the bearing axis. It is assumed in this model that the spring force is a 
linear function of the shaft stiffness (assumed constant) and magnitude of 
displacement. 

c. Damping Force - The damping force acts on the rotor center such 
that it resists translational motion by the rotor. The damping force is 
assumed to be a function of the damping coefficient experienced by the rotor 
in its particular medium and the linear component of the rotor's velocity. 

The nature of the rotor's dynamics readily lend the analysis to 
description using complex variables. Defining z as the time dependent 
position of the rotor's center, referenced such that z^ = 0 coincides with 
the rest condition where the rotor center is aligned with the bearing axis, 
then the equation of motion may be written as: 

£f = [-D z - K z ] = M z (1) 

L r c r c r g 



where z (t) = the vector position of the center of gravity of the rotor. 
S 

Furthermore, z (t) may be broken down into the vector sum: 



z(t)=z(t)+z (t) 

g c gc 



where z (t) represents the vector position of the rotor center of gravity 
gc 

relative to the disk center. 



The equation of motion then becomes: 
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or, 



M (z + z )=-Dz -Kz 
re gc re re 



Mz + Dz + Kz = Mz 
re re re r gc 



( 2 ) 



z (t) is simply a function of the radius of unbalance and rotational 
gc 

speed: 

N r jut 
z (t) = 6e 
gc 

Assuming synchronous whirl response, the motion of the rotor may be 
described by: 

. | | jut 

z c (t) " l z C l e 

Making these substitutions into Equation (2), and simplifying, results in 
the following: 



r u c"| 2 j(jjt 

LM 6 J co e 



z (t) = 1 

C(K p - u M p ) + j(uD p )] 

This harmonic motion of the rotor is termed forward synchronous precession 
since the orbital motion of the rotor is characterized by the forcing 
function's angular frequency and direction. Breaking down the rotor's 
mechanical impedance, 

[(K p - u 2 M p ) + j(uD r )] , 

into its magnitude and phase then leads to the harmonic response equation: 



, [M fi] w 2 e j(<0t "* ) 

z (t) = — 

c 2 2 2 ' ' ^ 

[(K p - u\) + (uD r ) ] 



( 3 ) 



The phase angle <j> describes the lag between the rotor's response and the 
forcing function. This phase lag which is due to damping is given by: 



